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In this paper we consider directed walks on a tree with a fixed branching ratio K at a finite 



temperature T. We consider the case where each site (or link) is assigned a random energy 

0' 

uncorrelated in time, but correlated in the transverse direction i.e. within the shell. In this 
paper we take the transverse distance to be the hierarchical ultrametric distance, but other 



possibilities are discussed. We compute the free energy for the case of quenched disorder 
and show that there is a fundamental difference between the case of short range spatial 
correlations of the disorder which behaves similarly to the non-correlated case considered 
previously by Derrida and Spohn and the case of long range correlations which has a totally 
different overlap distribution which approaches a single delta function about q = 1 for large 
L, where L is the length of the walk. In the latter case the free energy is not extensive in L 
for the intermediate and also relevant range of L values, although in the true thermodynamic 
limit extensivity is restored. We identify a crossover temperature which grows with L, and 
whenever T < T C (L) the system is always in the low temperature phase. Thus in the case of 
long-ranged correlation as opposed to the short-ranged case a phase transition is absent. 
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1 Introduction 



The problem of directed polymers in a random medium can be formulated on a lattice 
|H, H, ||, f|, [| , or in the continuum limit [§, [7|, |[ ||, |lOj, [TT|, [12] . On a lattice there is a random 
energy associated with each bond (or site). Walks (or polymers) start at a given point and 
are allowed to proceed only along the positive direction of one of the coordinates which is 
referred to as "time". The other coordinates are referred to as "transverse". The partition 
function is given by 



-PE(w) 



where the sum is over all walks w of L steps and 

E{w) = Yl € iJ 

(ij)Ew 



(2) 



is the sum of the random energies along the walk. (3 = 1/T is the the inverse temperature 
in the proper units. 

In the continuum limit the partition function is given by the functional integral 



Z{(5) 



J[Dx(t)]expl-P £ dt 



x(t) 2 + V(x(t),t) 



(3) 



where x(£) is the (d-l)-dimensional transverse position of the polymer at time t and V(x(t),t) 
is the random potential. The term (l/2)x(t) 2 measures the bending energy of the polymer. 
On the special lattice of a Cayley tree and with uncorrelated disorder 



(4) 



many properties of the model could be extracted analytically Jl3|, [14], like the free energy 
and the probability of overlaps between two walks. The model exhibits a phase transition 
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at finite temperature T c . Define q(w,w') to be the fraction of their length that two walks 
w, w' of length L spend together. The probability distribution for overlaps is then given by 

^) = (4EE%-<zK w/ )) exp(-(3E w - (3E W ,)\ . (5) 

For T > T c it was found that the probability of overlaps is a single delta function at q — 0, 
whereas for T < T c the probability distribution consists of a weighted sum of two delta 
functions: 

P(q) = ^K0) + (l - |) %-!)■ (6) 

This distribution implies that the free energy landscape consists of many valleys separated 
by large barriers. Two walks lying in the same valley have an overlap of q = 1, whereas 
walks lying in different valleys have zero overlap. 

In the continuum limit the model has been treated for general d by the variational 
approximation |J and by the l/cZ-expansion [II], both valid for large d. For the special 



case of d — 2 a Bethe ansatz technique yielded some exact results 0. In the continuum 
limit both the case of short-ranged and of long-ranged correlations of the disorder have been 
considered. The correlations have been defined by 

(V(x,t) V(x',t')) = 5(t-i'm(x-x') 2 ), (7) 

and classified as short-ranged or long-ranged according to the form of F(y 2 ). For the short- 
ranged case one usually takes || [15], [RJ 



F{y 2 ) = ^-{a, + yT^\ (8) 
7 - 1 

where g > is the strength of the disorder and 7 > 2 determines the range of the correlations. 
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The case of 1 < 7 < 2 is considered long-ranged. Also considered long-ranged are 
correlations of the form 

F(y 2 )=a - -^-y 2 ^\ (9) 
I-7 

with < 7 < 1. The constant ao is important to maintain the requirement \F(y 2 )\ < F(0) 
dictated by the Schwarz inequality for the appropriate y-range, but is sometimes neglected in 
the literature ||, likely because it contributes only a trivial constant to the free energy. In the 
case of a correlation of the form (||), a has to be very big, such that [oq(1 — 'y)/ gY ~ is 
greater than the system size. An example of this kind of correlations is the so-called random 
field case (see e.g. [|HJ), for which 7 = 1/2 and 

F(y 2 ) = h 2 (N-\y\), (10) 

where N is the system size. In the short-ranged case a phase transition in terms of the 
temperature (or the strength of the disorder) has been found, where the two phases differ 
in the nature of the distribution of overlaps. A one-step replica-symmetry-breaking solution 
has been found for T < T c when using the variational approximation. Both phases in this 
approximation were found to be characterized by the trivial wandering exponent v = 1/2 
defined by 



x 2 (L))ocL^ (11) 

(here the bar represents configurational or thermal average and the brackets refer to averaging 
over the disorder). This is unlike the exact analytical result at d — 2 which yields v — 2/3 
(superdiffusion) Jj], [J, and simulations that found z/>l/2forc?>2 |L7|, [Tj|. There were 



claims in the literature that the wandering exponent is greater than 1/2 for any finite d in 
the disordered dominated phase p^| . On the other hand there are claims that have gained 



more momentum recently [[19], [II, 2(J 2T, 2^], that the exponent becomes trivial (1/2) at an 
upper critical dimension, which is presumably d c — 5 (four transverse directions). 

For the case of long-ranged correlations the variational approximation for the continuum 
limit model predicts || 

(12) 



Flory 2(1+7) ' 

and it is not known if this result is exact at any finite dimension. In this case no phase 
transition is found as a function of the temperature (or the strength of the disorder) and there 
is an infinite-step replica-symmetry-breaking (RSB) solution (a la Parisi) for the appropriate 
order parameter, which manifests itself as a continuous (non-delta-function) part in the 
overlap distributions for the walks. 

The special case of long-ranged harmonic correlations (7 = 0) has been solved exactly 



by Parisi [53|. An inspection of his solution reveals that in that case the free energy is not 
extensive but grows as L 2 (where L is the size of the system in the temporal direction). 
Another exponent which is often referred to in the literature is the exponent characterizing 
the free energy fluctuations 

(F 2 ) - (F) 2 (xL 2 ". (13) 

The uj exponent is related to v through the scaling relation 

2 v = 1 + uj. (14) 

Thus for the harmonic case v = 3/2 and uj = 2. The harmonic case is special in the sense 
that its solution is replica symmetric. 

The case of long-ranged correlations, or in fact any non-zero correlations of the disorder, 
has not been investigated in the previous treatments of the model on the Cayley tree. The 



interesting results obtained for the long-ranged case in the continuum limit, and the fact that 
an exact solution has been found in the case of harmonic correlations, motivated us to carry 
out an investigation of the case of non-zero ranged spatial correlations on a Cayley tree, both 
for the case of short and long-ranged correlations. We obtain some interesting results which 
will be presented below together with some open questions. A thorough understanding of 
the directed polymer problem is particularly important because of its connection with the 
KPZ equation |M] and with the behavior of flux lines in high-T c superconductors |25|, ^B, 



27, . There is also a well-known mapping from the directed polymer problem to Burgers' 



turbulence, where the case of long-ranged correlations is of importance [29 



2 The tree problem 

We consider directed walks on a branch of a Cayley tree of coordination number K + 1 (see 
Fig. 1). Each bond branches into K new bonds in the forward "time" direction. For each 
site (or alternatively each bond ending at the given site) we choose a random energy e(t, z) 
where t designates the shell, < t < L, and z is a label within the shell that can take K l 
values. The random energies are chosen from a gaussian distribution satisfying 

(e(t,z)) = (15) 
(e(t,z)e(t',z>)} = 6 t ,« f(d(z,z>)). (16) 

Here f(d) is a function to be specified later, and d(z, z') is a distance among points (sites) 
belonging to the same shell. Energies at different shells are uncorrelated. 
In this paper we choose the following definition for d(z, z')\ 

d u (z, z') = number of steps for two walks starting at z and z' 

and moving backwards in time to meet. (17) 
6 



This is a hierarchical distance between two points. It also satisfies d u (z, z) = and < 



d u (z, z') < t within the t-shell. This distance is also referred to as an ultrametric distance [ 30 1 
(hence the subscript u), since it satisfies a stronger inequality than the ordinary triangular 
inequality, namely 

d u {z,z') < max(d u (z,y), d u (y,z')), (18) 

for any point y within the shell. This ultrametric distance is very different from an euclidean 
distance on lattices characterized by translational invariance in real space, but is sufficient for 
defining spatial correlations of the disorder within a shell, and make the problem amenable 
to an exact solution. 

A different choice for a distance which is more suitable for calculating the root-mean- 
square transverse distance can be defined as follows (TP]]. Let us label each branch of the 



tree by 1, • • • , K, which we call directions. These labels are a priori arbitrary on a tree, 
but once the choice is made it remains fixed at each branching point. For a given walk of 
length t starting at the origin we denote by Z\ the number of times the walk moves in the 
1-direction, by z 2 , the number of times it moves in the 2-direction, etc. We then associate a 
vector (z\—t/K,---, zk — t/K) with the end point of the walk on the t-shell. We denote this 
i^-dimensional vector by R(z). Note that there is not a one-to-one correspondence between 
points z on the tree and vectors R as there are different points which are associated with 
the same vector. The transverse distance between two points is defined as 

/ K \ V2 

d tr (z,z') = - i*CO) 2 J • (19) 

This distance always satisfies the inequality 

d tr (z,z') < V2d u (z,z'), (20) 
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for any two points z, z' . The advantage of this distance is that in the absence of disorder 
one has 

R = (21) 

w - iK 1 "*)*- (22) 

where the bar denotes configurational average over all walks of length L. Thus v = 1/2 as 
is expected for a random walk. In the presence of disorder with spatial correlations, this 
distance is harder to use in a calculation of the quenched free energy of the model, and 
further discussion of the use of this distance will be given in a future publication. 



3 The replica solution 

The method we use in this section is a generalization of the method used in Appendix 1 
of ref. \W\ for uncorrelated disorder. To calculate the free energy of the model we use the 



replica trick 



We can think of Z£ as the partition function of n different walks of length L emanating 
from the origin of a branch of a tree. Adopting Parisi's scheme for RSB in real space, we 
assume that the following arrangement of the n walks gives the leading contribution to the 
free energy in the n — > limit when L is large: 

• (a) The n-walks stay together for the first L(qi — q ) steps (where q = and < q\ < 
1)- 

• (b) The walks split into m\ bundles of (n/mi) walks each and remain so for a time 
L(Q2 - <?i)- 



• (c) Continuing in this way, in the j'th step the walks split into rrij groups each com- 
prising of (n/rrij) walks and remain so for a time L(qj + i — gj) = LAqj. 

• (d) Finally, at time tu = Lqu = L, the walks split into n individual walks. 
Thus 

= g < gi < ■ ■ ■ < Qm = 1 (24) 

1 = m < mi < • • • < Tom = n - (25) 

We also define 

x(qj) = n/rrij (26) 

and thus 

1 = x(q M ) < ■■■ < x(qj) ■■■< x(q ) = n. (27) 
For large L, (Z^) is given by 

M-l LA 

{ZD = maxmax [] (exp {-/3{n/ mj ) (e?> + • • -e% )})) % . (28) 

Xlji V 71 ]! j — Q 

Here ef denotes the energy encountered by a walker belonging to the i'th group at time t. 

The factorization in eq. (|28|) follows from the fact that random energies at different times 
are uncorrelated. The factor K mj is a geometrical degeneracy factor which is entropic in 
origin. It arises because each of the rrij groups can choose its next step among K different 
possibilities, and each possibility gives rise, as will become clear in the sequel, to a config- 
uration characterized by the same contribution to the partition function (same Boltzmann 
weight). One can also associate a combinatorial factor with different possibilities to assign 
individual walks to bundles when they split, but this turns out to give a total factor of n!, 
which becomes 1 in the n — * limit. 
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To proceed we use the fact that the random energies are chosen from a gaussian distri- 
bution satisfying equation (plj). It then follows that 

(exp{-(3(n/ mj )(e? + ---e%)}) = exp {^(n/m 3 -) 2 ((e? + ■ ■ ■ e% ) 2 )} , (29) 

((e? + ---eVf) = mj f(0) + 2j2N j/ f(£) (30) 

t 

where Nj^ is the number of pairs of distance £ at time t G among the rrij groups of 

walkers, where tj = Lqj. The coefficients Nj t £ satisfy 



m 



, + 2j2Nj, e = m). (31) 



A careful enumeration of the distribution of distances after t steps (tj < t < reveal 
that the following identity holds: 

N jMq .- qk)+A t = 7) (— 1 ^r) m l k = l,---,j (32) 

with At = t — tj. Since we are interested in the limit M — > oo, i.e. Aq — > 0, we will take 
iV^ to depend only on j and omit At in eq. (j3^) . In deriving eq. fl32|) we used the fact that 
the groups of walkers split at each time tj = Lqj according to the procedure described at the 
beginning of the section, as well as the definition of the hierarchical (ultrametric) distance. 
Substituting the result (|32| ) in eq. ( PD| ) we find 

(n/m,) 2 (( e W + ... e W) 2 ) = nx(qj)f(0) + n S>(g fc _i) - x(q k ))f(L(q 3 - q k )). (33) 

k=l 

In the limit n —>■ the inequalities (|25|) and (^7|) are inverted and hence 

= x(g ) < " " < <q 3 ) ■< x{q M ) = 1 (34) 



Using expression ([33|) in eq. (|29|) and subsequently in the formula for eq. fl2~8|) , we 

obtain 
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- \? £ ^ f ^ { Xi9k) A q k iqk ' l] ) f( L ^ - ( 35 ) 
where an extremum over x(qj) has to be taken. In the limit of large M, q becomes a con- 
tinuous variable in the interval [0,1], and x(q) becomes a function on that interval satisfying 
< x(q) < 1. The summations in eq. fl35|) become integrals and we have 

ZK = [ X ^L l nK + ±0* f dq x{q) f{0) ■ \fC dq f dp x'(p) f(L(q - p)), (36) 
L Jo x{q) 2 Jo 2 Jo Jo 

where x'(p) stands for dx/dp. This expression can be further simplified: First we perform the 
p integration by parts using the fact that x(0)=0. Next, in the second term which remains 
a double integral we change the order of integration 

/ dq dp = dp dq, (37) 
Jo Jo Jo Jp 

leading finally to a single integral. These steps yield 

ZEL = f 1 ^ i nK + \{? f l d qx{q) [f(0)-f(L(l-q))), (38) 
L Jo x[q) 2 Jo 

where x(q) is to be determined by extremizing this expression. Equation (^) is the main 
result of this section. 

Before we end this section let us mention a simple generalization of the problem. Equation 
(|16|) can be replaced by a more general form 

(e(t,z)e(t',z')) = 5 ttt , f(t,d(z,z')). (39) 

so the spatial correlation can vary with time (but energies at different times are still uncor- 
rected). The possibility of a time dependent width of the distribution of disorder in the 



uncorrelated case has been considered in ref. JT3|. Using eq. (|39|), we can repeat all the 
steps leading eq. (|36"D , and we find that it is replaced by 

=7^ = In If + \{? fdqxiq) f(Lq,0) 

L Jo x{q) 2 Jo 

-\P 2 fdq f q dpx\p) f{Lq,L(q-p)). (40) 
2 Jo Jo 
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It is easy to pursue this more general case by the same method presented in the next section, 
but we will not consider it further in this paper. 



4 Solution for the free energy and overlap distribution 



In order to extremize the expression for the free energy (^SJ) derived in the last section, we 
take a functional derivative with respect to x(q) to obtain 

In if 1 



x 2 (q) 2 



+ oW(0) - f(L(l - q)) = 0, (41) 



and hence 



Ty/2 In K . . 

x{q) = 42 

^//(0)-/(L(l-g)) 

This solution is valid for the range of q for which the inequality < x(q) < 1 holds. Otherwise 
one has to choose x(q) at its maximal (or minimal) allowed values. 

Before we consider some concrete candidates for the spatial correlation function, we 
comment about the posssible admissible correlations. The Schwarz inequality requires that 



the function f(y) defined in equation (jT6|) satisfy \ f(y)\ < 7(0). More generally for the k'th 
level the covariance matrix 

C z , z , = (e(k,z)e(k,z')) = f(d(z,z')), (43) 

must be positive semi-definite, which means that all its eigenvalues are non-negative. In the 
k'th shell the hierarchical distance takes values from 1 to k. Writing down the covariance 
matrix for the k'th shell we verified that a sufficient condition for an admissible correlation of 
the disorder is a function f(y) which is always non-negative and is a monotonically decreasing 
function of the hierarchical distace y. For simplicity we will limit the discussion for the case 
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K = 2, but it can be generalized to any value of K. First we present as an example the 
eigenvalues of the covariance matrix in the 3'd shell of the BL: 

/(0)-/(l), (multiplicity 4) 
/(0) + /(l)-2/(2), (multiplicity 2) 
/(0) + /(l) + 2/(2) -4/(3), 
/(0) + /(l) + 2/(2) + 4/(3), 

and these are all non-negative if /(0) > f(l) > /(2) > /(3) > 0. In general we now obtain 
by induction all the eigenvalues of the covariance matrix of the (k+l)'th shell once we know 
the eigenvalues of the matrix of the k'th shell. For k=l the covariant matrix is given by 



1 /(o) /(i) X 



(44) 



which has eigenvectors (1,-1) and (1,1) and corresponding eigenvalues /(0) — /(l) and /(0) + 
/(l). For k=2 the covariant matrix is a 4 by 4 matrix 



C(2) = 



(45) 



C(l) D(l)\ 
D(l) C(l) 

where C(l) is the k=l covariance matrix and D(l) is a 2 by 2 matrix, all the elements 



V 



/ 



of which are equal to /(2). The eigenvectors of C(2) are (1,-1,0,0) and (0,0,1,-1) with 
the corresponding eigenvalue /(0) — /(l) of multiplicity 2. In addition there are two new 
eigenvectors (1,1,-1,-1) and (1,1,1,1) with eigenvalues /(0) + /(l) - 2/(2) and /(0) + /(l) + 
2/(2). In the k'th level there are always k+1 distinct eigenvalues, the (k+l)'th of which is 
always given by 



Afc+i(fc) = /(0) + /(l) + 2/(2) + • • • + 2 k ~ l f(k) 



(46) 
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At level (k+1) the first k eigenvalues remain the same as level k first k eigenvalues, but with 
double multiplicity. This is because the corresponding eigenvectors are just obtained from 
the previous eigenvectors by adding zeros at the beginning or at the end. But there are two 
new eigenvectors (1, - - • , 1, — 1, • • - , — 1) and (1, • • • , 1) with corresponding eigenvalues which 
are given by 

X k+1 (k + l) = \ k+1 (k)-2 k f(k + l), 

h + 2(k + 1) = X k+1 (k) + 2 k f(k + 1). (47) 

It is thus straightforward to check that the positivity and (decreasing) monotonicity condi- 
tions mentioned above yield only non-negative eigenvalues. We now turn to some concrete 
examples: 

(i) The case of short-ranged spatial correlations. 
We consider first the function 



with A > 0. In this case eq. fl42|) becomes 



Ty2\nK/g 
K A -[a + L(l-g)]-A}i/2 



X ® = r-A_r „\ r,n_„M-A U /2 < 49 > 



and in the limit of large L the solution for x(q), < q < 1 becomes 

T <T r 



x(q) = { Tc (50) 



1 T>T, 



c 



together with x(0) = and x(l) = 1. Here T c is given by 



T c = JM-. (51) 
V21nir V ; 
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For T < T c , deviations of x(q) from the form given in expression fl5"0|) are of 0(L~ X ), except 
when 1 — q ~ 0(1/ L), and x(q c ) becomes equal to 1 for q c = 1 — 0(1/ L). 

The solution found above is the same as the solution found in ref. fl3| , for zero-ranged 
correlations. Thus we see that the case of short-ranged correlations is characterized by the 
same overlap function and free energy as the zero-ranged case. The overlap distribution 
function is given in terms of the solution for x(q) by |30 



m = (52) 

and hence 

m= \ + T<T^ (53) 

5(q) T>T C 
as discussed in the introduction. The free energy is obtained by substituting the expression 
for x(q) in eq. and we obtain: 



; j _ T c \nK + f(0)/(2T c ) = ^2f(0)\nK T <T C 
TlnK + f(0)/(2T) T > T c 

where we have dropped corrections on the r.h.s which vanish as L — > oo. Note that on 
the Bethe lattice with the hierarchical distance we did not find any substantial change as A 
crosses the value 1 as has been found in the continuum limit. 

(ii) The case of long-ranged correlations 

In this case we choose for the function governing the disorder correlation 

f(y) = a - g y a , (55) 

with a > and g > 0. We take a to be very large, so that it satisfies 



T « ^ < 56 > 
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for the range of temperatures we are interested in. Furtheremore, if L is such that it satisfies 

7) ■ (57) 

then f(y) is positive for the entire range of allowed distances and since it is a monotonically 
decreasing function of y it is a bona-fide correlation. The region of physical interest is a® — > 00 
with L large but still satisfying the condition Eventually though, if one is interested in 
the true thermodynamic limit L — > 00 one has to consider the case L > (ao/g) 1//Q . In that 
case we must define 

( a -gy a 0<y<(a / g y/ a 
f(y) = (58) 

[0 (a / g y/ a <y<L 

for f(y) to be a proper correlation function. We will comment about the thermodynamic 
limit later in the section. 

For now, starting with f(y) given by eq. (|55D with the condition ( |57|) being satisfied, the 
solution for x(q) (see eq. fl42|)) becomes 

x(g) = V^^T (59) 



? c < g < 1, 



with 



1 (T-2\nK\ l,a , v 

We can identify an L-dependent crossover temperature (at which q c = 0) given by 

tj - l) = yy" V2 (ei > 

For T fixed, as L becomes large, T C (L) grows with L and thus we find ourselves always in 
the low temperature phase of the model, characterized by q c > 0. (If on the other hand 
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T > T C (L) then x(q) = 1 for any q > 0). Assuming that L is large enough so T << T C (L) 
(but the condition (|57D still satisfied), we find by substituting the expression for x(q) in the 
formula for the free energy, eq. (j38|) 

_ f .aa"^,,.^^,,,.^ (62) 

which can be simplified to give 



- F = V2^ L i +Q/2 , (63) 
l + a/2 K ' 

where we dropped constant terms. We see that the free energy is not extensive for this range 
of L values, but rather proportional to L 1+Q / 2 . It is also temperature independent as it is in 
the low temperature phase of the short-ranged case. 

Let us now consider the distribution of overlaps. Using eq. ( |52D we find 



P{q) = < mr) 6 ^ + f ^wzfcz* - q - qc ; (64) 

q c < q < 1 

in the limit of large L this expression becomes simply 

P(q) = 5(l-q). (65) 

This becomes obvious by going back to the expression for x(q) which in the limit of large L 
become x(q) = for < q < 1 and a;(l) = 1, which can also be expressed as 



(66) 



Thus in the limit of large L the solution becomes replica symmetric for any fixed temperature. 

Before we close this section we show that in the true thermodynamic limit the free energy 
becomes extensive. To achieve this we must allow L > (ao/g) 1 ^ and use the correlation 

17 




function denned by eg. (|58|) . We still demand that a be very large so T « Jao/(2lnK) is 
always satisfied. In that case we find for x(q): 



where 



x(q) 



T/T c (a ) < q < q cX 

q c i<q< q C 2 



l <? C 2 < q < l 



T c {a ) 



q c \ = 1 



a 



2 In if 



dc2 = l- I 



1 {T 2 2\nK* 



1/2 



1/a 



1/a 



9 



In the limit L — ► oo one finds 



(67) 



(68) 
(69) 

(70) 



x{q)=T/T c (a ), < q < 1, 



(71) 



together with x(0) = and x(l) = 1. In the limit T << T c (ao) we see that practically 
x(q) = for < q < 1 and x(l) = 1 which amounts to P{q) = 6(1 — q) as before. The free 
energy though is given by 



F/L = J2a lnK 



(72) 



and it is thus an extensive function of L as it should be in the thermodynamic limit. 



5 Summary and discussion 



In this paper we have considered the case of directed polymers on a Cayley tree in the presence 
of correlated disorder. We have used the ultrametric hierarchical distance to introduce 



distance within each shell, and this distance is simple enough to enable us to solve the model 
exactly under the assumption of a hierarchical Parisi-type solution. 

We have found two different types of behavior depending on the range of the disorder cor- 
relations. In the case of short range correlations the solution behaves like the non-correlated 
case: There is a phase transition at a finite temperature and the two phases differ by the 
temperature dependence of the free energy and by the expression for the overlap distribu- 
tion which is non trivial in the low temperature phase. In the case of long range correlations 
there is no phase transition as a function of temperature (strictly speaking the transition 
temperature T c (a Q ) can be made as large as we please by choosing a to be large enough). 
This is similar to the behavior that has been found in the continuum limit ||. However 
we have identified an L-dependent crossover temperature which plays a role for a finite-size 
system. We have also found that in the large L limit the solution becomes replica-symmetric 
but with the overlap distribution peaked at q = 1 at any temperature, which is the case 
in the short-ranged case only at T = 0. This result is reminiscent of Parisi's solution p3| 



for the case of harmonic correlations in the continuum limit where no RSB has been found. 
This is in contrast to results in the continuum limit for the non-harmonic case || where the 
variational approximation yielded an infinite-step RSB for the long-ranged case (but also no 
phase transition). Another feature we have found in the long-ranged case, which is similar 
to Parisi's result [J23|, is the non-extensivity of the free energy in terms of L over a large 



range of L-values. Eventually as L — > oo, the free energy become extensive. The exponent 
1 + a/2 may be related to the exponent u which characterizes the free energy fluctuations, 
see eq. fll3|) , but this can be established only after carrying out a calculation of the free 
energy fluctuations on the Cayley tree. 

We should emphasize, that because of the non-euclidean nature of the hierarchical dis- 
tance on a tree, we could not establish a relation between the exponents A in eq. ( fj5|) or 
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a in eq. (|55|) to the exponent 7 denned in eqs. (H) and (|9|) which characterizes the range 
of disorder correlations on ordinary lattices embedded in euclidean space. Related to this is 
the fact that the separation between short and long range correlations occurs for 7 = 2 for 
ordinary lattices (see eq.(|8|) and ref. ||) whereas we find the separation to occur at A = 
(or a = 0) for the tree problem. 

There are various possibilities to extend this work further. One is to consider the trans- 
verse distance defined in section and attempt to solve the model, including the behavior 
of the root-mean-square transverse distance characterized by an exponent v. 

Also, we have only considered the case of a gaussian distribution of the disorder. Other 
distributions are possible, some better suited for calculating 1/d corrections (like the ex- 
ponential distribution) One can attempt to obtain these corrections for the case of 
long-ranged correlations of the disorder. 
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A BCDEFGH 

Figure 1: A branch of a tree with branching ratio K=2. The ultrametric distance between 
points A and B is 1, between A and C (or D) is 2 and between A and E (or A and F, G, H) 
is 3. 
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